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Thermoelectrics are important in physics, engineering, and material science due to their useful 
applications and inherent theoretical difficulty. Recent experimental interest has shifted to strongly 
correlated materials, where the calculations become particularly difficult. Here we reexamine the 
framework for calculating the thermopower, inspired by ideas of Lord Kelvin from 1854. We find 
an approximate but concise expression, which we term as the Kelvin formula for the the Seebeck 
coefficient. According to this formula, the Seebeck coefficient is given as the particle number N 
derivative of the entropy S, at constant volume V and temperature T, ^Kelvin = — {§§/} VT - This 
formula is shown to be competitive compared to other approximations in various contexts including 
strongly correlated systems. We finally connect to a recent thermopower calculation for non-Abelian 
fractional quantum Hall states, where we point out that the Kelvin formula is exact. 

PACS numbers: 72.15.Jf, 71.27. +a, 73.43.Cd 



I. INTRODUCTION 

A complete understanding of thermoelectric effects 
is important in the physical sciences where wide rang- 
ing applications utilize materials with large thermoelec- 
tric power S (Seebeck coefficient). Thermoelectrics of 
strongly correlated materials are of fundamental interest 
since they present an important and challenging prob- 
lem. Recent experiments have revealed that some ma- 
terials, such as sodium cobalt oxide Na^OoC^ (NCO), 
possess unusually large thermopower—, due in part to 
strong electron interactions 2 -. Frustrated systems', like 
NCO, might produce further surprises in enhanced ther- 
mopower in some situations 2 ' 4 . In addition, emerging 
work^ from the fractional quantum Hall effect (FQHE) 
is revitalizing thermopower as a tool to investigative the 
topological non-Abelian quasiparticles^ thought to exist 
at filling factor 5/2l. 

Here we present the Kelvin formula for thermopower, 
Skeivin- This is a formula inspired by Lord Kelvin's ther- 
modynamic treatment of this variable in 1854 s . It is 
found by reconsidering the sequence of taking the ther- 
modynamic and uniform limits, and is a valuable approx- 
imation to the exact, but computationally intractable re- 
sult, obtained via Onsager and Kubo's treatments"' 

For strongly correlated systems, such as the t-J model, 
^Kelvin is found to possess an accuracy between the rather 
coarse Mott-Heikes formulation, and a better argued high 
frequency limit formulation due to Shastryi and studied 
in Refs. Q and [III For intermediate couplings, such as 
the Hubbard model, we argue that ^Kelvin provides one 
of the best available approximations, it is better than 
the high frequency limit. In certain dissipation-less situa- 
tions, such as the FQHE, Skelvin is exact, thereby provid- 
ing an elegant and simple derivation for the thermopower 
formula used in Ref. [5j (derived originally in Refs. Il2i .) 



Skeivin is obtained by completing Shastry's argument 11 
for the "absolute thermopower" , i.e. S of an isolated sys- 
tem. Kelvin originally studied^ this object using the then 
available techniques, later he and others emphasized rela- 
tive thermopower between two materials. Let us revert to 
the absolute thermopower as a starting point and imag- 
ine a long isolated cylinder of material of length L subject 
to a time dependent electric field — V$ and temperature 
gradient VT. — V<1> couples to the dipole moment and 
VT couples to the moment of the energy density (cf. 
Luttinger— ). These fields individually generate a dipole 
moment linear in the fields to lowest order, and the con- 
dition for the cancellation of the two contributions, i.e., 
the zero dipole moment (or zero current) condition, leads 
to the thermopower S for a finite system size L at finite 
frequencies ui as S(L,u)) — ^^{L,uj}. 



The thermodynamic limit, L — > oo, and the static 
limit, to — > 0, must both be taken, as known from On- 
sager— and other o 13 i 14 . Kubo's exact formulas obtain in 
the fast or transport limit, where L — > oo before w — >• 0. 
Taking the static limit u) — > before L — > oo leads to the 
slow , where Kelvin's approximate formula arises and is 
expressible solely in terms of equilibrium thermodynamic 
variables. 



We transcribe this discussion to a more convenient 
periodic system, by trading the length scale L for a 
wave vector q x = 2n/h and the L — > oo limit by the 
uniform limit q x — > 0. The slow limit corresponds to 
Ym\{q x — > 0, lo — > 0}, and the fast limit corresponds to 
lim {uj — Y 0, q x — > 0}. The thermopower measures the in- 
duced thermoelectric voltage due to a temperature gra- 
dient and, as such, a useful and general formula for ther- 
mopower is given by the ratio between the thermoelec- 



2 



trical and electrical conductivities^, 

Xp(«.).*(-9.)( w ) 



S(q x ,u) 



where 



X A ,bW = * / dte^- 0+ )*<[i(i),B(0)]) 



(1) 



(2) 



(n|A|m)(m|B|n) (3) 



is the susceptibility of any two operators A and B, where 
p, K = H — /j,N, and J x are the charge density, the 
(grand) Hamiltonian, and the charge current operator, 
respectively, at finite wave vectors; H, [i, and N are 
the Hamiltonian, the chemical potential, and the total 
number operator, respectively. The susceptibility writ- 
ten in Eq. ([3]) is the Lchmann representation (where 
p„ = exp(—(3e n )/Z is the probability of the quantum 
state \n) with energy e n and Z is the partition function 
and j3 = l/k^T with fee the Boltzmann constant) which 
we find useful below. With Eq. (JTJ) , we can take different 
limits and obtain various interesting formulas. 



II. THERMOPOWER FORMULAE 



called Mott-Heikes (MH) approximation introduced by 
Heikes (popularized by Mott^) is written Smh = (a*o — 
/j,(T))/q e T, where fio = fx(T = 0). Smh is obtained by 
rather drastically replacing the first part of Eq. (Q| by 
the zero temperature chemical potential (j,q to make the 
theory sensibly behaved as T — > 0. From thermodynam- 
ics, we know that = {§^)e,v, and, hence, Smh 
relates thcrmopowcr to the partial derivative of entropy 
S with particle number N, at a fixed energy E and vol- 



ume V. We see below that Si 



Kclv 



is similar, but with 



The Kubo formula 



more natural "held" variables, namely T and V. 

C. High frequency formula 



From Eq. {!]), we can make a high frequency approxi- 
mation, where uj 3> w c (w c representing all finite charac- 
teristic energy scales), leading to the object S*. The for- 
mal expression and evaluation for S* are discussed else- 
where^ and we only quote the results. We have argued 
that S* is the best possible approximation to the exact 
Kubo formula for strongly correlated systems^ such as 
the t-J model, since the high frequency limit respects 
the single occupancy constraint and is closer to the DC 
limit than initially expected. It is not specifically suited 
for Hubbard type models, since the high frequency limit 
assumes w 3> U , and cannot capture the physics of cor- 
relations effectivelyii. We will see that SKoivin steps into 
this breach, and provides a very useful alternative for 
Hubbard type models^. 



Taking the fast limit and using the continuity equa- 
tions to pass from densities to current operators, Eq. (TTJ) 
gives the exact Kubo result^ 



S 



Kubo — 



1 J °° alt ^ dr{jZ(t - it)J x {Q)) M (T) 

f °°dttfdT(j x (t- i T)j x (o)) q e T- 



T 



(4) 



where q e is the charge of the carriers and J E the energy 
current. 



B. The Mott-Heikes formula 

For narrow band systems, such as NCO, high Tq 
superconductors, or heavy fermion systems, the so- 



D. The Kelvin formula 



To obtain an approximate thermodynamical expres- 
sion, we consider the slow limit of Eq. (fTJ). S is among 
the few objects (along with Hall constant and Lorentz 
number) where this process gives finite and approximate 
results, unlike the electrical conductivity where the slow 
limit gives meaningless results^. 



This limit is identified with Kelvin since he essentially took the equilibrium limit of an interacting gas of particles. 
The slow limit (q x —> 0, lu — > 0) is easiest to compute starting from Eq. (TJJ 

Sk**n - Hm i'^-'-^i - (5) 



= urn : — ; : ■ — — — . 

9--^rXp(«.),p(-a.)(0) 



To simplify we first consider the numerator of Eq. (|S|) which we rewrite by first using the Lehmann representation and 
then taking the q x — > limit. Note that p(q x ) tends to a conserved quantity q e N and cannot mix states of different 
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energy so e m -)■ e n . Thus, 



lim > 



Pn - Pn 

e 



(n\p(q x )\m)(m\K(-q x )\n) 



m c n 



Pn - Pr, 



(n\p(q x )\m)(rn\k(-q x )\n) 



-m °n 



1 _ e P(e n -e m ) 

lim } p n {n\p(q x )\m)(m\K(-q x )\n} 



£m^£n Em En 

n.m 



= <le^2 pp n 8 en>em (n\N\m)(m\k\n) 



q e /3[(NK)-(N)(K)} 
d , » , d , - , 



(6) 



The derivative with respect to p in Eq. ((5]) is within the grand canonical ensemble and performed with a fixed V 
and T. The denominator of Eq. (J5j> is treated similarly yielding q 2 (3[(N 2 ) — (N) 2 } = q 2 d(N) j dp. Combining it with 
Eq. ©, yields 



1 il(H)-p£(N) 



lvin 



dfi 



(7) 



To further simplify Eq. (|7|) we note a relation found 
in textbooks on thermodynamics in the grand canonical 
ensemble: (H) = E = fl + TS + pN (CI the grand po- 
tential), so that (ff ) T ,v = m($^)t,v + r(g) T ,v and 
hence 



^ Kelvin 



g e \dNJ TV 

-1 /v 

9e \ 9T / A . 1 



(8) 
(9) 



where we used, to go from the second equality to the last 
equality, a Maxwell relation obtained with dF = —SdT — 
pdV + pdN, and equating g TdN = g Ndx - We refer to 
the last two equivalent equations (Eqs. (J5]) and ([5])) as 
the Kelvin formula for the thermopower. This formula is 
unknown in the literature as far as we are aware. 

Note that Smh is similar to ^Kelvin- The distinction is 
that in ^Kelvin, the number derivative of the entropy is 
taken at constant T rather than at constant E. Thus, in 
the low T limit of a metal, where p{T) oc T 2 , they differ 
in the linear T coefficient by a significant factor of two. 
We show below that for non-interacting electrons, scat- 
tered by impurities, ^Kelvin is closer to the exact result 
than Smh- Further, we see that the approximation of ex- 
changing the slow and fast limits has some justification 
in dissipation-less systems, such as in the FQHE where 
Skeivin i s identical to that found by several workers (see 
below) . 



III. APPLICATIONS OF THERMOPOWER 
FORMULAE 

A. Free electrons 

To gain insight into the strengths and weaknesses of 
the various thermopower formulations discussed above 
we consider non-interacting degenerate electrons treated 
within the limit of elastic scattering at the Born level 
with an energy momentum dependent relaxation time 
t(p,uj). This is a modestly dissipative system, but at 
such a simple level that the Boltzmann-Bloch equation is 
an adequate description. The solution for S is available 
in textbooks and a useful benchmark for various approx- 
imations. In the low temperature limit^i, to 0(T 3 ), 

ir 2 k 2 d 

SMott = T ~^ M "jJ L ln [P^){(Vp?T{p, p)) a ] \a^a . 

(10) 

a formula often ascribed to Mott and po(p) is the sin- 
gle particle density of states per unit volume per spin. 
In this non-interacting electron context, S^eivm gives (to 
0(T 3 )), 



S K 



clvin 



= T- 



w 2 k 2 d 



In [po(p)} 



fill 



3<? e dp 

which differs from the exact answer (Eq. (JTOJ) ) in the 
neglect of the relaxation time r and particle velocity 
in the logarithm. Smhj to the same order, gives 



TT 2 k 2 d 

Smh = T——^- — In [p (p)} \a^u 
6q e dp 



(12) 
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which is off by an important factor of two from Sxcivin 
(Eq. (fTT]0. The formulations (Mott-Heikes and Kelvin) 
would be identical if // oc T, which occurs if the system 
possesses a ground state degeneracy, and in the classical 
regime. The high frequency approximation gives a better 
result than all these and, again in the low temperature 
limit, to 0(T 3 ), 

s* = T ^^ ln h^)(K) 2 U I™ ■ (13) 

Other than the neglect of the energy derivative of t, this 
is the same as the exact result. Hence, ranking the ther- 
mopower approximations for non-interacting electrons we 
have, from worst to best, Smh, SKeivin, and S* with the 
exact result being Smou- 



B. Hubbard model 

For intermediate coupling models, the relative rankings 
of the various approximations can be different. In partic- 
ular, ^Kelvin can be superior to S*, since the effect of cor- 
relations is diluted in the latter by making the assump- 
tion of to ;§> U , whereas S^civm retains ui <C U. The sign 
of the true (i.e. transport) thcrmopower and the trans- 
port Hall constant are expected to flip as we approach 
half filling in the Hubbard or t-J models due to the on- 
set of correlations (carriers become holes measured from 
half filling rather than from a completely filled band) . In 
the case of the t-J model, the high frequency Hall con- 
stant R* H and S* do display this behavior—. However, 
for the Hubbard model, R* H and S* do not display a sign 
change 16,17 . SKeivin °n the other hand, does appear to 
show the expected change in sig n 16 ' 18 . Further discussion 
concerning the relative merits of SKeivin and S* will be 
reported later—. 



C. NCO and the t-J model 

To show the usefulness of S'Koivin, we apply it to NCO 
since (i) we have previously investigated- this system 
while benchmarking S* , (ii) the system is intrinsically 
interesting 1 , and (iii) we can compare different ther- 
mopower formulations on equal footing. As discussed 2 , 
the action in NCO takes place primarily in the cobalt ox- 
ide planes where d-shell spin-1/2 electrons live on a trian- 
gular lattice and these strongly interacting 2D electrons 
can be modeled with the t-J model. Hence, we exactly di- 
agonalize the t-J model on a L = 12 site two-dimensional 
triangular lattice with periodic boundary conditions (cf. 
Fig.[U). Note that we only show results for the t-J model 
with zero super-exchange interaction (J = 0), as the re- 
sults only weakly depend on J. To map the t-J model 
to NCO we follow Refs. i and! and give results as a 
function of electron doping x = |1 — n\ away from half 
filling (n is electron number density). 




i i i i i i i i i i i i i i i i i i i i i i 

2 4 6 8 100 2 4 6 8 10 



771*1 771*1 

FIG. 1: Thermopower vs. T for the t-J model (with J = 0) 
corresponding to NCO in the Curie-Weiss metallic phase near 
x ~ 0.7. (a) and (c) correspond to x = 0.67 and x — 0.75 
for the NCO system (i > 0) while in (b) and (d) the sign 
of the hopping has been switched to investigate the enhance- 
ment expected for frustrated systems. The black, red (light 
gray), and blue (dark gray) lines are 5*, Skdvm, and Smh- 
Finite size effects at low T are treated in the spirit described 
previously^. At each x, for T below an appropriately cho- 
sen cutoff temperature To = 0.5|t|, the thermopower is fit to 
S(T) -> aT + bT 2 where a and 6 are obtained from the com- 
puted S(To) and S'(To) providing a sensible extrapolation to 
low T and plotted as dashed lines. The inset figure (e) depicts 
the 12-site unit cell. 



S* adequately describes the physics of NCO for x > 
0.5 and, in particular the so-called Curie- Weiss metallic 
phase 2 near x ~ 0.7. The subject of this work, however, 
is Skeivin- We see in Figs. QJi and c and [2^,, similar 
to SmHj Skeivin does a good job capturing the physics 
with minimal computational effort. However, Skoivin 
does seem to overestimate the thermopower for inter- 
mediate temperatures and high dopings as compared to 
Smh- Near x ~ 0.7, Skeivin and Smh are similar but as 
x is decreased the two formulae diverge and for low dop- 
ings, Skeivin better captures the physics as it is closer to 
the more accurate high frequency limit S*. 

An interesting property of the triangular lattice un- 
derlying the physics of NCO is its geometrical frustra- 
tion^, cf. inset Fig. [Tfe. It was predicted 2 -^ that if the 
sign of the hopping amplitude were flipped to t < the 
thermopower would be enhanced at low to intermediate 
T. We have considered this situation in Figs. [TJd and d 
and (2b. Since the thermopower enhancement for t < 
compared to t > is largely a consequence of electron- 
electron interaction it is important to determine whether 
this effect is captured by SKeivin- We see this enhance- 
ment is captured to some extent by Sxoivin and Sxeivin 
is better than Smh in the large doping region where the 
enhancement is the greatest, but is missing some of the 
electron-electron physics at very low T that is captured 
by S* (as is Smh)- 
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FIG. 2: The thermopower vs. T and x for the t-J model 
(with J = 0) for t > (a) and t < (b). Note that x ~ 0.7 
corresponds to the Curie- Weiss metallic phase of NCO, cf. 
Fig.[T] The line- type and color-coding is the same as in Fig.[T] 
Note that for dopings below 0.5 it is not clear whether the t-J 
model adequately describes the physics of NCO. 

D. FQHE at v = 5/2 

We now discuss how ^Kelvin i s applied to dissipation- 
less systems such as the FQHE where thermopower can 
be used as a possible non-Abelian quasiparticle detector—. 

For a weakly disordered electron system (from Eq. (|10p 
and (jlip) Skdvin essentially gives the dissipation- less 
thermopower where particle velocities are further approx- 
imated. If the system is dissipation-less and the parti- 
cle velocities are also energy independent, such as the 
FQHE, then we expect Skeivin is exact. An expression 
for the thermopower in a 2D electron system in the pres- 
ence of a perpendicular magnetic field (the FQHE sys- 
tem) has been derived 5 -^, assuming zero impurities, as 
(Eq. (6) in Ref. Hj). Yang and Halperin show 5 - that 

S ~ ksN \og{d) where d > 1 is the quantum dimension 



of the quasiparticles for the FQHE at v = 5/2 (provided 
they are non-Abelian). Thus, a non-zero entropy linear 
in N is obtained. From Eq. (JHJ, we see that the ther- 
mopower is the derivative of the entropy with respect 
to the number of particles at constant T and V . When 
entropy is linear in particle number, as in non-Abelian 
FQHE states, dS/dN — > S/N and the formulas are iden- 
tical. Our derivation provides a simple and straightfor- 
ward insight into the formula given previously 5 . 



E. High Temperature Superconductors 

Before concluding, we point out an intriguing applica- 
tion of S'Kcivin for high T c superconductors. For different 
families of high T c compounds, a universal curve of the 
thermopower, at T = 290 K , as a function of hole concen- 
tration p ~ 1 — n has been observed^. The thermopower, 
in all families, vanishes near optimal doping [p ~ 0.16) 
starting out positive at small p. Phillips et al. 20 appeal 
to the atomic limit of Smh as an explanation. Viewing 
this data — more generally, through the prism of S^civm 
(Eq. ((HJ) we conclude that the optimal filling, i.e., maxi- 
mum T c , additionally corresponds to a local maximum of 
the electronic entropy as a function of filling. This con- 
clusion is powerful, since we avoided the difficult issue of 
calculating either thermopower or entropy, merely using 
the link between them provided by ^Kelvin- 



IV. CONCLUSION 

It is clear that Smb., which has served as a virtual 
workhorse for years, has a new competitor in S^eivm- 
This simple minded approximation can be written in 
closed form and in many difficult regimes, where the ex- 
act Kubo-Onsager expressions are not useful, and Skeivin 
provides an excellent guide to the physics of the system. 
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DE-FG02-06ER46319 at UCSC, and MRP acknowledges 
support from LPS-CMTC at University of Maryland and 
Microsoft Project Q. 
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